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Abstract 

In this paper, we study the bead model: beads are threaded on a set of wires on the plane 
represented by paraUel straight hnes. We add the constraint that between two consecutive 
beads on a wire, there must be exactly one bead on each neighboring wire. We construct a 
one-parameter family of Gibbs measures on the bead configurations that are uniform in a certain 
sense. When endowed with one of these measures, this model is shown to be a determinantal 
point process, whose marginal on each wire is the sine process (given by eigenvalues of large 
hermitian random matrices). We prove then that this process appears as a limit of any dimer 
model on a planar bipartite graph when some weights degenerate. 



1 Introduction and presentation of the bead model 

We consider the collection of configurations of beads strung on an infinite collection of parallel 
threads lying on the plane. A bead configuration on these threads gives a configuration of points 
on Z X R. We impose the following constraints on the configurations: 
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Figure 1: A bead configuration. 



• The configuration must be locally finite : the number of beads in each finite interval of a 
thread must be finite. 
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• Between two consecutive beads on a thread, there must be exactly one bead on each neigh- 
boring thread. 

Let fl be the set of bead configurations satisfying these two conditions. 

The first aim of this paper is to construct probabiUty measures for our infinite system Q that 
are uniform in a certain sense. 

If there were only a finite number of threads of finite length, and a fixed number of beads on 
each thread, then the set fl would be a bounded convex set of R^, where is the total number of 
beads. Therefore the normalized Lebesgue measure on would give a uniform probability measure. 
Therefore we look for probability measures that satisfy the two following properties: 

• they are ergodic under the action of Z x M by translation 

• conditioned in an annular region, they induce the uniform measure on allowed configurations 
inside this region. 

Such a probability measure is called an ergodic Gibbs measure. When endowed with an ergodic 
Gibbs measure, the set fl is called a bead model. 

This model can be viewed as a interface model on I?. Indeed, the position of the beads have the 
same combinatorics as the square lattice I? (rotated counterclockwise by 7r/4). A configuration can 
therefore encoded by a height function : — > M, where (f>{x, y) is the ordinate of the bead {x, y) 
on its thread. The problem of existence (and uniqueness) of Gibbs measures for this model can 
thus be formulated in terms of random surfaces with a simple attractive potential |l5[| reflecting the 
hard-core interaction between beads. However, this approach does not lead to an explicit expression 
for the Gibbs measures. We adopt another point of view that will allows us to give a closed formula 
for cylinder events. 

The (T-algebra of events for our probability measures is defined as follows. To each bounded 
Borel set i? of Z x M and to each bead configuration a> G f2 is associated an integer Xsioj). Let "J 
be the smallest a-algebra such that all the maps Xb : — > N are measurable. 'J is generated by 
the elementary events 

{cj e 17 I Xb{uj) = n]. 

If P is a Gibbs measure on (il,?"), it defines through the application X : B ^ Xb a random 
process with value in the set of boundedly finite, integer- valued measures, i.e. in other words, a 
random point field. 

Even without giving for the moment any explicit description for a Gibbs measure, it is possible 
to estimate, at least heuristically, the probability of rare events. For example, one can estimate the 
probability that n beads lie in the same wire interval of length e. For this event to occur, then, 
due to the geometrical constraint imposed on configurations, there must be n — 1 beads in a small 
interval of size e on the left neighbor thread, and n — 2 beads in the same interval on its left, and 
so on, and the same must happen on the the right side of the considered thread. That is why the 
probability of this event must be of order e"+2Efc=i ("-fc) = _ Note that this is much smaller 
than the probability of having n points of a Poisson process in such a small interval, which is of 
order e": there is thus a kind of repulsion between beads, induced by the geometrical constraint on 
bead configurations. 

Such a repulsion has been observed in some point processes on the real line, especially in 
determinantal point processes |p^ , for which correlation functions are expressed as determinants of 
a certain kernel. Certainly, the most famous example of such processes is the so-called sine process, 
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which describes the statistics of the eigenvalues of large random hermitian matrices with Gaussian 
entries (GUE ensemble in the bulk of the spectrum, and whose kernel is given by he following 
expression: 

sin(a; — y) 



k{x,y) 



it{x - y) 



We will see that the Gibbs measures we construct on ($7, J) define determinant random point 
fields on Z X R, for which correlations functions are given by determinants of a kernel J, whose 
restriction to a single thread is the sine kernel. 

We prove the following theorem: 

Theorem 1. For a fixed average density of beads, there exists a 1-parameter family of ergodic Gibbs 
measures (P^) on (f2,?'). When endowed with one of these measures, is a determinantal 

random point field on Z x R, and the probability of any cylinder event has a completely closed form. 
In particular, the marginal on each thread is the sine random point field. 

The exact expression for the kernel is given in Theorem 2. The parameter is directly related to 
the average distance between a bead and its right neighbor just below it and describes the amplitude 
of a magnetic field that tends to push the beads in some directions. 
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Figure 2: Two typical bead configurations for difli'erent values of 7. The parameter 7 is negative in panel (a), and 
positive in panel (b). 

A way to construct these Gibbs measures is first to consider a discretized version of the bead 
model. The set of possible configurations fit C 51, is constituted by the configurations for which 
the beads are located at sites of a lattice with mesh t. We show in sections 2 and 3 that in this 
discrete setting, there exist probability measures supported by fit, for which the distribution of 
the beads is a determinantal random point field, by exhibiting a bijection between the discretized 
bead model and the dimer model on the honeycomb lattice — or equivalently random tilings of the 
plane by rhombi. The measures on random tilings have the Gibbs property. Then we prove that 
the sequence of discrete determinantal processes indexed by t converges to a determinantal random 
point field on Z x M when t goes to zero. In the second part of the paper (section 5 and 6), we prove 
that the bead model appears not only as a limit of the dimer model on the honeycomb lattice, but 
as the universal limit behavior of the dimer model on any bipartite periodic planar graph. 

The fact that the correlations between dimers is given by the (discrete) sine kernel has been 
noticed before in particular cases by Johansson |^ and Baik et al Our approach gives a 
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systematic method to find a subset of edges, for wliicli correlations can be made arbitrarily close 
to the sine process. 

2 A discrete version of the bead model and the dimer model 
on the honeycomb lattice 

We assume for the moment that the threads are not continuous lines, but a one-dimensional lattice 
with mesh size t. The possible positions of the beads are labeled by coordinates 

{x,ty) eZxtl, (1) 

X representing the thread on which the bead lies, and y being the coordinate running along the 
thread. 

There is a mapping between discrete bead configurations and dimer configurations on the hon- 
eycomb lattice H. Before explaining this mapping, we recall some facts about the dimer model 
on H. 



2.1 The dimer model on H 

A dimer configuration of the honeycomb lattice H is a subset of edges of H such that every vertex is 
incident with exactly one edge of this subset. A dimer configuration is also called in graph theory a 
perfect matching. On the set of all possible dimer configurations, there is a two-parameter family |^ 
of Gihhs probability measures, all elements of which satisfy the following properties: 

• They are ergodic under the action of the lattice translations Z x Z, 

• If a dimer configuration is fixed in an annular region of H , the dimer configurations on the 
inside and in the outside region are independent, and the dimer configurations inside are 
uniformly distributed. 

Here is how these measures are defined, following ^ . The vertices of H are colored in white 
and black such that no two neighbors have the same color. Weights are assigned to edges of H 
according to their orientation: a, 6, c. A fundamental domain of H is obtained by taking for 
example a white and a black vertex sharing an horizontal edge. The vertices of H are labeled by 
their color and the coordinates of their fundamental domain {x,y). 

The local statistics for the Gibbs measure corresponding to these weights have a determinantal 
form: the probability that edges ei — (wxi^yi^h^i^ yi^), . . . ^e^ = (Wa-^. j,^., bj,/^ j,^) belong to the 
random dimer configuration are given by 

Hb,;.,,,w,^,,J] (2) 

where K is the so-called Kasteleyn operator. In the case of the honeycomb lattice, it is simply 
the weighted adjacency matrix restricted to the lines corresponding to white vertices, and columns 
corresponding to black vertices: if w and b are neighbors, then K(w, b) = a, 6 or c depending on 
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the orientation of the edge (w, b). If they are not neighbors, then K(w,b) = 0. Since K is 1? 
periodic, its inverse can be expressed using inverse Fourier transform over the unit torus: 



K ^{hx',yl,'Wx,y) = K ^(b(^/_2,)^(j,/„y), Wq^o) 



T2 a + b/w + cz/w 217: z 2iiTw 



(3) 



When none of the weights is greater than the sum of the other two, one can show (by computing 
expHcitly the integral above) that every type of edge appears in the random dimer configuration 
with positive probabihty. The measure is then said to be liquid |^ . 

For such weights, there is particularly well-adapted embedding for the hexagonal lattice, the 
so-called isoradial embedding Q corresponding to these weights, defined as follows: the length of a 
dual edge of H is equal to the weight of the corresponding primal edge, so that the dual faces of H 
are represented as triangles with side length a, b and c. The dual faces for the usual embedding of 
H with regular hexagons are equilateral triangles and the this embedding corresponds to equally 
distributed weights a = b = c. 

The measure and the embedding (up to a global scaling factor) depend only on the ratios a/b 
and c/b. For our purpose, we will chose a = t, & = 1, and c = exp(7t). These values satisfy 
the triangular inequalities for 7 e (—1,1) at least for t small enough. The distance between two 
successive horizontal edges is t. We denote by P^.t and K;^ ^ the probability measure and the inverse 
Kasteleyn operator corresponding to these particular weights. 

2.2 Correspondence between beads and dimers 

The mapping we construct between discrete bead configurations and dimer configurations can be 
described as follows: there is a bead at (x, ty) if and only if the horizontal edge incident with the 
white vertex in fundamental domain [x, y) is in the dimer configuration. 

A way to see geometrically this correspondence is using these isoradial embeddings of the hon- 
eycomb lattice described above. Take an isoradial embedding of the honeycomb lattice for weights 
a = t, 6 = 1, c = e'^'^ for 7 G (~1, 1)- Once chosen a dimer configuration on iJ, draw a bead in the 
middle of an horizontal edge if it appears in the dimer configuration and you end with a discrete 
bead configuration with mesh size t. Reciprocally, from a bead configuration one can reconstruct a 
dimer configuration by placing horizontal dimers on edges crossing an occupied site, and completing 
the configuration. This is always possible, because of the intertwining of bead positions. Moreover, 
the completion is unique as soon as there is at least one bead on each wire. 

For a fixed t, each value of 7 corresponds to a liquid Gibbs probability measure on dimer 
configurations, that can be transported to bead configurations. The local statistics of the beads 
coincide with that of horizontal dimers. The probability measure on bead configurations benefits 
from the conditional uniform property of the dimer Gibbs measure. 

This procedure defines for a given t a family parameterized by 7 of probability measures on 
discrete bead configurations that have the conditioned uniform property. The correlations between 
beads are given by determinants: the probability of having a bead at the sites (xi, tj/i), . . . , (xfc, ty^) 
in the random bead configuration uj is 



P^,t[(a;i,t?/i),...,(xfe,tyfc) e w] =r det K ;(b, 




(4) 
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where ^ is defined by 



K-(b.„w)=// 



3 Construction of explicit Gibbs measures for the continu- 
ous bead model 

In this section, we wih give an exphcit description for Gibbs measures for the bead modeL But 
before proving theorem 0, it is necessary to investigate the behavior of the kernel defining the 
discrete bead model. In other words, one has to compute the asymptotics of K:^ t(b^,w) when t 
becomes small and y large. 

The first thing to note is that for the weights we chose, the probability of an horizontal edge at 
a given white vertex (that is the probability of a bead at a given site in the discrete model) is 

P,,, [horizontal edge] = tK-^dWo, wo,o) - t J J^^ . t^T^^ + o{t). 

(6) 

In order to keep a constant average density of beads, we must choose the rescaled vertical coordinate 
5 equal to tyy^l — 7^. 

The asymptotics of the kernel for this vertical scaling are given by the following lemma. 

Lemma 1. In the vertical scaling limit t — + 0, tyy/l — 7^ ^, the coefficients -^-^=L=K^^^{hx,y, w) 
converge to 




e-^'^'f'i-t + ifpV'^-l^fTr- ifx<0. 

zvr 



1,1] 
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1 f - ^i'^(^) 



In particular, when x — 0, 

Proof. As we said before, the entries of the inverse Kasteleyn operator are given by 

/, Nil z'^w^ dz dw 



2 t + - (1 + zeT*) 2i7rz 2i7rw 



To evaluate this integral, we first perform the integration over w by the method of residues. If 
a; > 0, the rational fraction 

" 7 — rrn — ^ 

tw + (1 + zeT^^) 

has one pole a,t w = wq{z) = 

By Cauchy's theorem, the integral 



77- [ fz{w)dw 



is zero unless the pole 100(2) is in the unit disc, i.e. 

Re{z) < -i±^;l^ = -1 + (1 - j')t^ + 0{t'). (9) 

Define 6*0 = 6*0(7, t) = Arccos (^^^^^r^) = t^I -"f^ + 0{t'^). The constraint (|) on the pole 
to be inside the unit disk can be rewritten as 

arg(z) e (vr - 6'o,7r + 6*0) 



Posing z = — e*^ = _g»t0v^^)^ j-^ integral, we get 



K- (b.,„w) = / z-y f-l±Ii:lY ^ (10) 



'Ro(z)< 

e-'.»(^^l (11) 



2eT 

^e'^'^e*'' - 1\" d6l 



V / / ~ity<h I 



27r 



- 1 yr^d^. (12) 



In the vertical scaling limit t — s- 0, ty^l — 7^ — > ^, we have 



lim = 1, lim e-'^y^"/^- ^ e"^'*' , lim ''''^^^"^ ^ ^ = 7 + i0 V^" 



tyr^7^ 



7" 
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Thus, the integral above, muhiphed by ^. , converges to 

1 



K-(b,,„w) = — / e-'«^ (7 + ^071^)" 




When a: < 0, has two poles: there is a pole at w = in addition to that located at 

w = Wo (z) = -i±f^. Since wf,{w) goes to zero when z goes to infinity, the sum of the residues 
is zero. Therefore the integral of fz{w) on the unit circle is not zero only if wo{z) is outside of the 
unit disc. It equals in that case the opposite of the residue at ^0(2). Again with the change of 
variable z = — e*^ = _g«t0V^-7^^ -^yg have 

K.-Mb..,w)^-/ (13) 

(14) 

= I / + j 1 e-.^V-.^ I ejy^^V^^\ ' y/T^d^^ ^^^^ 

/l-'y2 t\/l--y2 y \ / 

Thus in the scahng hmit, by Lebesgue dominated convergence theorem, 

lim^^lE=K-(b,^„w) = -i / e-'«^(7 + *0)"d0 

what terminates the proof of the lemma. □ 

From the exact expressions ( p^ and (|l5| ) of K"^(bif w), one can easily check that the entries 
are uniformly bounded in y and t for a given value of x, leading to the following lemma: 

Lemma 2. V a; G Z 3M^>Q y t < to VyeM \K-\{h^^y,w)\ < M^. 

These two lemmas will now be used to prove theorem 0, stating that this family converges 
weakly to the determinantal random point field on Z x M with kernel J-y. 

Theorem 2. For each value ofj G (—1, 1), the discrete bead model converges weakly when t goes to 
to a determinantal random point field on Z x R. The kernel of this limiting determinant random 
point field is J-y : 




1.11 ^ ^ 



— ifx>0 
2tt 



e-'^"^ {-i + i(t>y/l--iA — ifx<0. 
[-1,1] ^ / 27r 



The marginal of the process on a given line is a determinantal random point field on M with kernel 
It is thus the sine random point field of the eigenvalues of large random Hermitian matrices. 
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Proof. Since tightness is automatic for random point fields Q], it is sufficient to prove the conver- 
gence of finite dimensional distributions in order to prove the weak convergence of the family of 
random point fields (fi, 3^, P^^t)- 

Let Ii, . . . ,Ik be segments on wire xi, . . . ,Xk respectively. It will be convenient to use multi- 
index notations 



Wrijl, \n\ 



3 = i 



We will prove the convergence of the moment generating function G^j^^{zi, . . . , Zk) of the joint 
law of (X/j ; . . . , Xi^ ) 



11(1 --.)^^^- 



n 



(18) 



(19) 



where E^.t is the expectation with respect to the probability P-^.t on discrete bead configurations. 
The factorial moments 

are quite easy to compute. They are given by the formula: 

Aill{ni,...,nk) = E P^,t [there are beads at (a;i,ty^ ),..., (ccfcjty^ J] (20) 

y}...y^ G — /-*- distinct 



where the sum is performed over all the distinct n^tuples of /j, j = 1, . . . , /c. By equation (Q), this 
can be rewritten in terms of determinants of matrices with blocks of size ni, . . . ,nk 











tl"l det 








distinct 









■ (21) 



V^-'-vi, g / , distinct 



l<«l>il<ni 



The determinant in the sum is the same as the one obtained by multiplying which converges when 
t goes to zero to 



4^)(n)= / det 





















det 








d"e 





















(22) 



where the integration variable S, is the n-tuple , ■ ■ ■ , Cm i ■ ■ ■ : ^nj)- Since the coefficients of ^ 
are bounded uniformly in t and y, say by M, then using Hadamard inequality, we get a uniform 
bound on the coefficients Ail^^(ni, ■ ■ ■ ,nk) 



\A[%)\<Y[\i,r^{^\M) 



|n| 



(23) 
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Therefore, by an argument of dominated convergence, the entire series (3^j,^t(2),z € C'"' converges 
uniformly on compact sets towards 

OW(z)= ^4^)(n)i-^ (24) 

which is the moment generating function for the hmit distribution of {Xj-^ , • ■ • , Xi,,). The probabihty 
of having for ail j G {1, . . . , k} exactly Uj beads in Ij is given by the following formula 

('—I'll"! f)"- 

P, [X,, = n„ ...,X,,=n,] = L-i__gW(z)|^^(^^^^^^^^. (25) 

In particular, the probability of having no bead in a Borel set B is given by the Fredholm 
determinant 

[XB=0]=Bct(ld-XBK-\B)=Q^{l)^y2^—^ det[j^(6-0)]d"e. (26) 

„=o -^s" 

□ 



Many quantities about the continuous bead model can be easily computed, using the underlying 
dimer model. An example of such a quantity is the average ratio between the distance between a 
bead and its neighbor on the left and above it, and the distance between two successive beads on 
the same thread. This average ratio is just the limit of the proportion of c-edges among the non 
horizontal edges in the random dimer configuration of H. It is then equal to 



V r 11 7,t [c-edge] e'^'^K \,(-l,-l) arccos7 

M - 1- P,,, or 6-edge] = I'^o K-^-l, 0) + e^^K-^-l, -1) = ^''^ 



This quantity would have been difficult to obtain directly from the description of the Gibbs measure 
for the point process, but has a simple interpretation in terms of dimers. 

4 Comments and interpretations of the bead model 

4.1 GUE matrices and uniformly distributed intertwined points 

Afterwards, it may seem not so surprising that these configurations of "uniformly" intertwined 
points are related to the determinantal sine process and more generally to random matrices from 
the GUE ensemble. 

Take a random (finite) matrix Hn from the GUE ensemble, and define for k g {l,...,n}, 
Hii^ — {Hn)i<i,j<k to be the submatrix of formed by the first k rows and k columns of Hn- 

(k) (k) ~ (k) 

Let A]^ < Aj, be the eigenvalues of Hn ■ It follows directly from the mini-max formulation for 
the eigenvalues that 

V fc e {1, . . . , n}, V J- e {1, . . . fc} Af < Af ) < Xf^\'\ (28) 

Moreover, a result of Baryshnikov Q states that if we condition on the eigenvalues Ai, . . . , A„ of 
Hn, then the eigenvalues of the submatrices are uniformly distributed in the simplex 

§(Ai,...,A„) = < [Xj )l<k<n~l \\j <X) ' < Xj + i ] X) ' < X) < X)^^ ^ . (29) 
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The bead model is somehow a bi-infinite analogue of this. 



4.2 The bead model as an asymmetric exclusion process 

A bead configuration can be interpreted as the history of a collection of particles located on a one 
dimensional lattice Z and jumping from left to right. Time is continuous and is flowing vertically 
along the threads, and there is a lattice site between two successive threads. Joining every bead to 
the bead just above it on the neighboring right thread, one gets an infinite collection of monotonous 
paths representing the trajectories of the particles: a bead on a thread corresponds to a jump of 
a particle from the site at the left of the thread to the site on its right. Because of the geometric 
constraint on beads, these paths cannot touch each other. Consequently, the particles are submitted 
to an exclusion rule: a particle cannot jump to a site if this one is already occupied by another 
particle. 



Figure 4: The trajectories corresponding to the bead configuration of figure ^ and positions of particles (black 
squares) at different times ti and t2. 



The Gibbs measures on bead configurations, viewed as families of monotonous paths con- 
structed as above, are probability measures on all possible evolutions of particles. The Gibbs 
property and ergodicity imply that the marginal of these measures for a fixed time ( i. e. along an 
horizontal line) give stationary measures for some Markovian dynamics. 

The discrete bead model gives a discrete version of this particle system: in this picture a particle 
is represented by a c-edge, and a hole by a 6-edge. Under P-y, the average particle density p is equal 
to the limit of the probability of a c-edge, and therefore related to 7 by the following expression 

p = lim Pt.^ [c-edge] = 1 - !^££^ (30) 

SO that the density is an increasing function of 7. 

The Asymmetric Simple Exclusion Process (ASEP) is also an example of particle systems with 
the same constraint of exclusion. Its evolution is Markovian, and the transition rates from an 
allowed configuration to another is constant. The translation invariant stationary measures for this 
model are Bernoulli probability measure, whose parameter is the density. 

If the particles are not located on the vertices of the finite lattice Z but on a finite annulus 
Z/7VZ, then the number of particles is a conserved quantity. For a fixed number of particles, the 
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stationary measure is uniform for ASEP. This is not the case for the bead model with a finite 
number of beads"'^ the probabihty of a configuration of particles depends not only on the number of 
particles, but also on their positions. 

The properties of the particle system coming from the bead model differ from that of ASEP. It 
would be interesting to study more in details these properties using the dimer microscopic structure, 
and to compare them with that of ASEP, that are also related to random matrix theory fl^. 



5 The bead model as a universal limit for dimer models 

Although the bead model was presented in the last section as the limit of the dimer model on the 
honeycomb lattice, it turns out to be much more general. Indeed, the bead model appears as the 
limit of any dimer model on a planar Z^-periodic bipartite graph. We first recall some facts from 
the theory of the dimer model on a planar bipartite lattice(see P, 0| for more details). 



5.1 The dimer model on a bipartite planar periodic graph 

Let G be a planar bipartite Z^-periodic graph, together with a periodic weight function on the 
edges of G. We suppose that the fundamental domain, delimited by an horizontal path and a 
vertical path 7j,, contains n black vertices bi, . . . , b„ and n white vertices Wi, . . . , w„, and that G 
has at least one dimer configuration. 

There is a two-parameter family of Gibbs measures on dimer configurations of G for these 
weights, parameterized by the two component of an external magnetic field {Bx,By) One can 
associate to this weighted graph, as in the case of the honeycomb lattice, a Kasteleyn operator K 
that will describe the dimer model on G. For a given value of the magnetic field B, the probability 
that some edges ei = (wi, bi), . . . , = (wfe,bfc) appear in the random dimer configuration is 
given by the following formula 

[ei, . . . ,efc] = I n K(w„b,) I det [K-(b„w,)] (31) 

where the entries of K^^ are given by the inverse Fourier transform 

KB^(bU,w^)= // [K{e^^z,e^yw)]-'^^= [[ Q^AfUz^A^^ (32) 
JJr^ 2i7Tz2iTrw JJj2 P{eB-z,e^yw) 2iTrz2iTrw ^ ' 

The characteristic polynomial P(z,w) is the determinant of the Fourier transform K{z,w) of 
the periodic operator K, and Q{z, w) is the comatrix of K{z, w). The asymptotics of K^j^ and thus 
the correlations decay depend on the regularity of the integrand in (^) , and in particular on the 
presence of zeros of P{e^''z, e^^w) on the unit torus. 

The spectral curve { (z , w) S | P{z, w) = 0} is a complex algebraic curve of a special kind : it 
is a Harnack curve ^ ^ For generic values of {B^, By), P{e^^z, e^^w) has zero or two roots 

bead model for a finite number of threads can be constructed following the same procedure as in the beginning 
of this article. We impose the number of threads TV to be even to ensure that the geometric constraint on beads 
makes sense. We get leading to a 1-parameter family of determinantal random fields, whose kernel is obtained by 
replacing the integral in (fTj) by a discrete sum. 
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on the unit torus, and the phase diagram describing the behavior of the measures in function of 
and By is given by the amceba of the spectral curve, i.e. the image of P{z, = by the mapping 

Log : {C*f 

{z,w) ^ (loglz^oglu-l). 

When {Bx,By) lies in the interior of the amoeba, the characteristic polynomial P{e^'' z,e^^w) 
has two conjugate roots on the unit torus and the correlations decay polynomially, and the corre- 
sponding measure is said to be liquid. When {Bx,By) lies inside a bounded connected component of 
the amoeba, the correlations decay exponentially fast, and the measure is gaseous. When [Bx,By) 
lies in a unbounded complementary component of the amoeba, the measure is said to be solid and 
there are infinite deterministic dual paths crossing no dimers with probability 1. 

The existence of a dimer configuration on G is equivalent to that of a dimer configuration on 
the torus Gi = G/1? . We suppose that such a configuration on Gi exists, that can be lifted to a 
periodic dimer configuration Cq of G. A dimer configuration C can be interpreted as white-to-black 
unit flow, i.e. a 1-form with divergence +1 at each white vertex, and divergence -1 at each black 
vertex. Therefore, the difference C — Co is a divergence-free flow. The horizontal (resp. vertical) 
slope of a Gibbs measure ^ is the expected amount for n of the flow C — Co across ~jy (resp. 72;). 
Two Gibbs measures having the same slope are in fact equal ^. The Newton polygon N{P) 
of P, the convex hull of the exponents of monomials of P, coincide with the set of all possible 
slopes for a Gibbs measure on dimer configurations. The structure of this amceba is related to the 
geometry of N{P) through the Ronkine function 

H : iB...B,) „ II log(P,e«...e«.»,) Ji-^. (33) 

In particular, the solid and gaseous phases are mapped to integer points of N(P). 

To clarify all these notions, we apply them to the particular example of the honeycomb lattice we 
discussed before: choosing weights a, b, c for edges according to their orientation (without magnetic 
field) is in fact equivalent to choosing all weights equal to a = t and imposing a magnetic field equal 
to Bx = log(c/6) = 7t, By = log(a/6) = logt. The fundamental domain of the honeycomb lattice 
is constituted by one black vertex and one white vertex. Therefore, the Fourier transform K{z, w) 
is 1 X 1 matrix:if (z, w) =t(l-t-l/w-|-z/u;). The characteristic polynomial P{z, w) is therefore also 
equal to t(l + 1/w + z/w), such that 

P(e^^z,e^^w) + -(1 + ze'^''). (34) 
w 

The cofactor Q{z, w) of K{z, w) is by convention equal to 1. The Newton polygon and the amoeba 
for this model are represented on figure ^. 

The bead model is obtained from the dimer model on the honeycomb lattice when the magnetic 
field goes inside one of the outgrowths of the amoeba, in the thin region separating two solid phases. 
In the general dimer model, the bead model will also appear near the frontier between the liquid 
phase and solid phases. But before explaining how to find the bead model in this setting, we need 
some more information about the local geometry of the amceba, in particular about its unbounded 
outgrowths of the amceba: the tentacles. 
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(a) 



(b) 



Figure 5: The dimer model on the honeycomb lattice; (a) the fundamental domain of the graph, (b) the Newton 
polygon N{P) of the characteristic polynomial P{z,w) = 1 + (1 + z)/w, and (c) the amoeba of P in the plane 
{Bx,By). 

5.2 Tentacles of the amoeba 

Consider a particular side of the Newton polygon N{P). Changing the generators of the 7? lattice 
acting on G by translation induces a linear transformation of N{P). A change of basis of 7? is 
encoded by an element M of SL2(Z). The linear transformation acting on N{P) is {M~-^)'^. After 
possibly such an operation, we can assume that this side is horizontal, and that the polygon lies 
above it. Recall that the Newton polygon represents the possible slopes for a Gibbs measure on 
the dimer model on G. When the slope of a Gibbs measure is a lattice point of the boundary of 
N{P), the system is a solid phase. 

We want to investigate the geometry of the phase diagram for values of the magnetic field 
inducing measures with a slope close to the particular side of N(P) we chose. In particular, we 
seek for the shape of the curve in the neighborhood of the frontier between the liquid phase and 
the different solid phases, corresponding to the points of the particular side of the polygon. 

To get a measure with a slope close to that side of N{P), we apply to the system a magnetic 
field oriented essentially downward (Bx,By) = (c, —R). To remain close to the notations used in 
the previous section, we introduce the small parameter t = e~^. 

When t is small, the leading terms in the characteristic polynomial P(e'^z,tw) are those with 
the smallest power in w, say Sq. 



By a suitable choice of the origin of the Newton polygon (deforming the paths and/or jy 
delimiting the fundamental domain of G), one can assume that Sq = and that all the roots of 
■fo(^) = Y.-yO'foX'' are positive ||]. 

If e'^ is not a root of ^o(-^) = J2-y (^■yoX'^ , then for t small enough, P{e'^z, e^^w) has no roots 
on the unit torus. In this case, the magnetic field {Bx,By) = (c, — i?) belongs to a unbounded 
component of the amoeba. The corresponding measure fi{Bx,By) is frozen. On the contrary, if 




(35) 
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e'^ is a root of this polynomial, then for every R large enough, the polynomial has two complex 
conjugated roots on the unit torus: we are in the liquid phase. The amoeba defining the liquid 
phase has therefore tentacles going to infinity with asymptotes given by the straight lines x = c. 

For generic weights, the asymptotes are all distinct, and there is one asymptote for each segment 
between to lattice points on the side of N{P). Moreover, one can give an asymptotic expansion for 
the equation of the boundary of the amoeba: since the boundary of the amoeba is the image of the 
real locus of the curve, it is given by the equation 

P{±e^-,±e^y) = (36) 

In the neighborhood of {Bx,By) — (c, — oo), admits an asymptotic expansion in t = e^^": 
Bx ^ c + at + O(t^). Since Pie", 0) = 0, we have 

P(e^^±e^«) =P(e=+^i"+°(*'\±t) = (cle=alP(e^ 0) ± 92P(e^ 0))t + 0(t2) (37) 

Therefore the coefficient ci in the expansion is defined by 

_ d2P{e^,0) 
''-^e^d,P{e^,0) 

and the two curves B^ = c± ^~'' dlp[e'''ff) define asymptotic branches of the boundary of the 
amoeba in the neighborhood of (c, —oo). Define (3 as 

_ d,P{e^,0) 

9lP(e^0) ^^^^ 

For any 7 e (—1, 1), the curve 

P. = c + 7/3e^« (40) 
lies inside the amoeba for By negative enough. 



5.3 Deep inside a tentacle 

5.3.1 Analytic results about the roots of P 

Let us fix c to be equal, as before, to the logarithm of one of the roots of Pq- For a fixed z, the 
polynomial P{e'^'^'*^^ z, W) has d roots Wo(z), . . . , Wd^\[z). Since e"^ is a root of Pq, one of these 
Wj(z\ say it is Wq, equals when z — \. If all the roots of Pq are distinct, Wq(z) is the only zero 
having this property. The 2-to-l property shows that Wq(^z) does not equal zero for z ^\ 
Therefore, by compactness of there exists an e > such that 

Vje{i,...,d-1}, Vze§i W,{z)\>t (41) 

Differentiability of the roots with respect to the coefficients and symmetry by complex conju- 
gation imply that there exists ^(t) = QqX + O(t^), with 6*0 > such that 

\Wq{z)\ < t ^ z e [e-'^(*),e'^(*)] (42) 

In fact, an expansion of P similar to ( ^7| ) shows that when arg(z) ~ 0(t) 

M^o(^) = 7t + ^^^ + 0(t2), (43) 



and therefore, 6l(t) t/3^1 - 7^ + 0{\?). 



15 



Fi gure 6; A tentacle (in thick lines) and curves described above, for different values of 7. 



5.3.2 Asymptotics of the inverse Kasteleyn operator inside a tentacle 

The Newton polygon N{P) can also be obtained indirectly from the max flow-min cut theorem as 
the intersections of half-planes made of points satisfying linear constraints to be a slope of measure 



The side of the Newton polygon we are looking at is a segment of one of the straight lines 
delimiting one of these half-planes. When the average slope of the Gibbs measure lies on this line, 
one of the inequalities defining N{P) becomes an equality, implying the existence of frozen dual 
paths, with a direction perpendicular to the side of N{P), that are not crossed by any dimer with 
probability 1. These possibly frozen paths will be the threads of our bead model. When the slope 
is not exactly on that boundary of N{P), some dimers may cross these paths. We will see that 
these defects will play the role of beads strung along these threads. 

Let e = (w, b) be an edge of G crossing one of these threads. When the slope of the Gibbs mea- 
sure is on the side of the Newton polygon, this edge appears in the random dimer configuration with 
probability 0. In particular, there is no dimer configuration on the torus Gi = Gjl? corresponding 
to a lattice point of the side of N{P) containing this edge. As the cofactor Qe{z, w) = Qhw{z, w) is 
up to a sign the determinant of the Kasteleyn operator on Gi \ {w, b}, it contains only monomials 
of degree at least 1 in w (otherwise, there would have been a Gibbs measure with a slope on the 
side of N{P) for which e is a dimer with positive probability, what is in contradiction with the fact 
that e crosses a frozen path). 

We determine now the asymptotic expression for the coupling function K;:^\(bx, w) correspond- 
ing to our magnetic field {B^, By) ^ (c + /37t, logt), between w and the black end bx of a translate 
ex of e by X = (x, y) € J? 



on dimer coverings 




(44) 
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Proposition 1. Denote by J^(x,^) the kernel of the bead model. 




e-*«"^(7 + i0Vl -72)^^ ifx>0 
M^,0 = {-''-Y' .^^ ^ !"d0 (45) 

/r\[-i,i] '^'^ 



In the scaling limit t — > 0, , t(3yy^l — 7^ ^, the coefficients t(bx,w) have the following 
asymptotics 

K-Ubx,w) ^ tpeVl-7'^a(x,0, (46) 
where the quantity pe is given by 

Once multiplied by K^, it represents the proportion of this type of edges among the defects along a 
thread. 

Proof. We denote by f{z,w) the rational fraction inside tlic integral ( |4^ ) 

, z-y-^w^-^Q4e^+'''>^z,tw) 
= P(e-+^7t,,t^) ■ 

The integral defining \ is evaluated by performing first the integral over w. Suppose first 
that X > 0. There is no singularity at w = 0, since deg^Q > 1. The only pole in the unit disc 
is Wo{z)/t when z e [e~*^^*^ e*^^*^] = I^. In this case, we introduce ({(f>) = 7 + i(f>\/l — 7^ with 
z = e'l^'^l'V^^ and get: 



K~ubx,w)=y^ 



z-y{Wf,iz)/tf 'Q,{e-+l'-''z,Wo{z)) dz 

d2P{z,Woiz))t 2iTTZ 



l+0(t) g-»/3ty0yi3^(^(^) ^ Q^(gC+t/3C(0)^^^(^)^o(^2))^^^ 



^ 7-i+o(t) a2P(e^+*'3C(0),tC(0)+O(t2)) 27r- ^ ^ 

For a small t and a fixed we have 

g^(gC+t/3C(0)^ tC(<^) + 0(t2)) = tC(0)a2Qe + ©(t^), (50) 
52P(e=+t/3C(0)^ tC(0) + 0(t2)) = + 0(t). (51) 

where the derivatives of polynomials P and Qe without specified point are evaluated at (6*^,0). For 
the first expansion, we used the fact that Qe{z, 0) = and therefore diQ^ — 0. Therefore, 
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When a; < 0, the rational fraction in the integral has a multiple pole ai w = which is hard to 
evaluate directly. However the rational fraction is a o{^) when \w\ — > oo, and hence the sum of all 
residues in the plane is 0. 

Let us bound the residues at the simple roots of P: Wi{z), . . . ,Wd-i{z) . We know already 
from (^ ) that there exists e such that for every j G {1, . . . ,d—l}, \Wj{z)\ > e for every z e S"'^. By 
the same argument of compactness, there exists a constant M > such that for t small enough, 



Vj e {!,..., d- 1} VzeS^ 



\d2Pie^+^^^z,Wjiz))\> and 



< M (52) 



and therefore, 



{Wj{z)/trQ,{e^+^f^\Wjiz)) 



< (e/t)^Ar = 0(t-^) 



Wj{z)d2P{e''+'^l^^z, Wj{z)) 
Thus the contribution of these residues is negligible as soon as x < —2. In that case, we have: 

f dz f dz 

K-^^(h^,w) = J^^Res^^a f(z,w)— + Res^=wo(.)/t f{z,w) — 

= -V'/ ^es^=w,(.)/,fiz,w)^~[ Res, 



(53) 



dz 

=w^M/^f{z,w)— (54) 



Using the estimates ( pO| ) and (|5l|), we find, using the same change of variable z — e*'3*viT^'^ 
as above, that 



Res„=wo(.)/t/(z, w) — 



z~y{^^Y-'Q,[e-+-<P'z, Wo{z)) dz 

2iTTZ 



2^7^ 7si\/, td2P{e-'+il^^z,Wo{z)) 

I -l + o(l) U3|t\/l- 



e~'^^"'^\/^(C(0) + Q(t))^-i 
t92P(e=+'5^«("*),tC(0) + O(t2))'^ 

Qe(e^+*^«(^\tC(0)+O(t2))g. (55) 



By Lebesgue's dominated convergence theorem, and by|50| and 51, the integral is asymptotic in 
the scaling limit t ^ 0, tPy^/l — ^ ^ to the following expression 



d2P 



and thus, 



K:;:t(bx,w) 



e-*^«(7 + z<^yT3^)-^ + 0(1) 



(56) 



(57) 
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When X = —1, the residues at the poles Wi{z), . . . , Wd-i{z) are not neghgible any more. How- 
ever, in this case, the pole at w = is simple. A direct evaluation of the integral shows: 



f dz f dz 

K:^,t(b(-i,y),w) = y ^ Res,„=o/(2:,'i«)^ + J ^es„=wo(.)/^fiz,w) — 



2i7r 

2 , 



2-^t92Qe(e=+''^*^,0) dz f z-y{^^^) ^Q,{e^+^^^z,Wo{z)) dz 



p(gc+/37t;j^ 0) 2iTTZ Jj^ td2P{e''+f^'^''Z,Wo{z)) 2iTTZ 

Posing again z = e''^*'^Vi and C(0) = 7 + i<p\/l — 7^, one has 



(58) 



K-Ubx,w)-t2/?yT^ J 



i+o(i) g-i/3v'tyVi^Q3(e'=+/3<(«,0) d0 



-'-i+od) Wo(e*'''^Vi-^')2a2P(ec+/3tC(0)^O) 
Using the following estimates from Taylor's formula 

d2Q{e'+'^^^^^\Q) = asQe + 0(t) (60) 

= - ta<P){d2P + 0(t)) = -tC(0)92P + 0(t2) (61) 

together with (|5^) and (^l|), and applying Lebesgue dominated convergence theorem after an in- 
tegration by parts, one can prove that in the scaling limit, the coefficient of the inverse Kasteleyn 
operator Kz^^^ih^^w) is asymptotic to 

ftO f , , , ,^ n 7^,-1 d(^ 



t/^yr^^ / e-'*«(7+i<^yr^) 



^2^' yR\[_i4] 277 



□ 



Remark. The ratio d2Qe/d2P controls the density of the limiting bead model. If one plugs the 
value (/) = 1 into ( pO[ ) and (|5l]), one can see that it is, up to terms of higher order in t, equal to 

*Qe(e^+'3^^zo,t«io) ^g2) 



id2P{e''+l^f^zo, tWQ)tWQ 

where (zo,wo) are zeros of the characteristic polynomial P{e''^'^''^-,t-) on the unit torus. When 
multiplied by Kg, the numerator is the length of the dual edge e* in the natural application from 
the dual graph G* to R'^ described below while the denominator is that of the vertical side of the 
fundamental domain. 

Lemma 3 (^). Let (e^^'zo, e^^wo) be a root of the characteristic polynomial, with {zq,wq) on the 
unit torus. The 1-form 

e = (w,b) ^ iK^h{e^- zo,e^ywo)Qh^{e^'' zo,e^ywo) (63) 
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is a divergence-free flow. Its dual is therefore the gradient of a mapping from G* ~ C. 

This mapping 5** is Z'^ -periodic and the symmetries of its range are generated by 

X = ie^^ ZQdiP{e^^ zq, c^^Wq) and y = ie^ywod2P{e^''' zq, c^^wq). 

This application ^> coincide with the notion of isoradial embedding for dimer models with 
critical weights (see 0) and gives a geometry well-adapted for the study of liquid measures on 
dimer configurations (see for example [Q). 

Proposition |l| shows that the kernel giving the correlations has the same form as the original 
bead model. However, in order to recover the bead model, one can not just look at one type of 
edges on threads but at all of them. 

Since the frozen paths have been chosen to cross no dimer when the slope is on the side of the 
particular Newton polygon we are looking at, they are bordered by white vertices on their left, and 
black vertices on their right. For a reason of parity between white and black vertices, there is no 
dimer configuration of the graph Gi deprived of the projection of these two vertices having a height 
change^ on the side of N{P). Therefore the arguments of the proof of proposition ^can be applied 
to obtain similar asymptotics as those given in that proposition for the coefficients of ^^^t between 
these vertices. 

Proposition 2. Let h^^y and w be respectively a black and a white vertex each bordering one of 
these paths, and in fundamental domains separated by a lattice translation {x,y). In the scaling 
limit 

t^O, t/3yVl-7'^e, (64) 
the coefficient K^^;(b2:_j^, w) has the following asymptotics 

^^^t(^x,y,^) ^ tp],^J^{x,£,), (65) 

where 

a2Qbw(e^0) 

These coefficients pi-,w are in fact the product of two terms, one depending only on b and the 
other on w. This property is stated in the following lemma 

Lemma 4. The rank of the matrix d2Q{e'^,0), restricted to projections of vertices bordering a 
thread, is equal to 1. In particular, for any h and any w bordering a thread, there exist Uh and 
such that 

Pbw - C/bKv. (67) 
Proof. The matrix d2Q{e'^,0) is the limit when t goes to zero of 

-^Q{e'+^'''zo,two) (68) 
two 



^The term height change here is an abuse of notations, since the difference between the reference unit flow Co 
and the one corresponding to any dimer configuration of Gi deprived of the two vertices has a non-zero divergence. 
What we mean here by this expression is in fact the powers in z and w in the weight of the configuration computed 
using the magnetically altered Kasteleyn operator K{z, w) divided by that of the reference dimer configuration. 
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where (zq, wq) are zeros on the unit torus of P^e'^'^^'^^z, tw). These zeros depend on t and 7 and 
their the first term in their expansion in t is obtained by plugging cj) — ±1 into equations ( pO[ ) and 




zo = l + 0(t), u.o = 7±Vl-7' + 0(t). (69) 

Since (e'^^'^'''*zo, two) is not real, it is a simple zero of P = det _fir, since the mapping from the 
spectral curve to its amoeba is 2-to-l out of the real locus. Q is the comatrix of K, its rank is 1 
at a simple zero of P. Therefore, as the limit of sequence of rank 1 matrices, the matrix 92(3(e^, 0) 
has a rank a most 1. As there is at least an non-zero entry in this matrix, the rank is equal to 1. 

The coefficient pbw is a multiple of the entry (b, w) of this matrix. Its decomposition into a 
product comes from the representation of a rank-1 matrix as a tensor product of a vector and a 
linear form. □ 



5.4 Convergence to the bead model 

We already said that the threads of our bead model would be the infinite collection of vertical paths, 
translated one from another, that are frozen, i.e they do not cross any dimer when the slope of the 
measure lies on the boundary of the Newton polygon. The beads are represented by the dimers 
crossing these paths when the magnetic field lies in one of the vertical tentacles of the amoeba. 

Like in section ^, as the magnetic field goes deeper into the tentacle of the amoeba, the picture 
of the graph in the plane is rescaled in such a way that although the probability of seeing a dimer 
crossing these "almost-frozen" paths goes to zero, the average number of such edges by centimeter 
of thread stays almost constant. The scaling limit we perform is 



0, 



To find the limiting distribution of this beads, we first evaluate the quantities 



(70) 



(71) 



where Xj. is the number of dimers crossing the (rescaled) thread interval. We look in details at 
the case fc = 1 when only one thread interval is at stake. The other cases are similar. For a given 
n, and a fixed value of 7 and of the scaling parameter t , we have 



E, 



7:t 



Xi\ 



(Xi - n) 



E 



■ 7.t 



[ei, . . . ,e„ e e] 



(72) 



ei,...,e„G-r 
distinct 



where the sum is over all possible n-tuples of edges crossing the thread interval /. The edges 
crossing / are labeled by their type (i.e. their projection on Gi = G/l?) and the coordinates of the 
fundamental domain they belong to. The edge represents the edge of type j in the fundamental 
domain with coordinates x = {x,y). The type label j ranges from 1 to d. Since the probability of 
having two such edges in the same fundamental domain is negligible, we can rewrite this sum of 
probabilities, as a sum over the fundamental domains and the types of edges crossing the thread 
interval. 

Xil 



E 



7,t 



(Xi 



- nV 



E E 



- 7,t 



,eiiee]+o{t) 



(73) 



xi,...,x„e/ jl 
distinct 
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The different probabilities P^_t • ■ • ,6^'^^ e C] are given by the determinant of a n x n matrix 
that is equal to 

{tp^l^T ( n X [p^.M^k - xutp^l^{yk - yi))] + 0(t"+i) (74) 

\l=l / - ' ~ 

Since is the product of two terms UjVk, by carrying out of the determinant by n-linearity these 
coefficients, equation (Q) becomes 

(t/3yT^)« [fl kA X det [U,,V,,.Mxk - xi,tp^/T^iyk - yi))] + 0(t"+i) = 

\l=l / - ' - 



(t/3v/l-7^)" n ^^'U^'^i X , det [J^ixk - x,,t/3v/l-72(yfe - yi))] + 0(t"+^) = 

\i=l / ^ 

(tp^l^r ( n ) X det [Mxk ~ x,, t/3yT^(yfc - y{))\ + 0(t"+i) (75) 



V/=l 

Summing now over the different types of edges, one gets 

d 

J2 P^,t[ei\,...,ei;ee] = 



KjPj is the proportion of edges of type j along the thread. These coefficients sum up to 1, and we 
have finally that expression ([73|) is a Riemann sum for 



/■■■/„ i<t:kn ["^^ ^"^^ -^i^ik-^i)] d" c (77) 

and the same argument of domination as in the proof of theorem || implies that the distribution of 
Xj converges to the distribution of beads in the interval / in the bead model of parameter 7. The 
generalization to any finite dimensional distribution is notationally cumbersome, but straightfor- 
ward. These considerations give thus the proof of the following theorem 



Theorem 3. Let 7 G (—1,1). In the scaling limit t — )■ 0, tf3yy^l — 7^ the point process 

describing the position of rare edges on the threads, identified with the almost frozen paths converges 
to the bead model of index 7, i.e. the determinantal point process on Z x K with kernel J-y. 

Recall that 7 describes the different possible ways to go deep into a tentacle. This theorem 
states that the bead model, with its 1-parameter family of Gibbs measures, is the universal limiting 
behavior of any dimer model on a bipartite periodic planar graph when the order parameters 
{Bx, By) go to infinity in staying in the liquid phase. 
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(a) 



(b) 



Figure 7: (a) A 3 x 3 fundamental domain of the honeycomb lattice. Weights of edges around white vertex w^j 
are a^j, bij and Cij. The Fourier multipliers have been distributed over the edges by gauge transformation, (b) The 
amoeba for generic weights on the graph represented on the left panel. This amoeba presents a gaseous phase, and 
three tentacles in three directions, each corresponding to a collection of frozen paths drawn on figure H. 



6 Interaction between bead models 

It often happens that a side of the Newton polygon is not the result of a unique frozen path, but 
that different paths give the same constraint on the slope. In that case, we do not have just one 
family of frozen paths, but several parallel families of thread, carrying all in the scaling limit a bead 
model. In this section, we describe the interaction between these different bead models in the case 
of the generic case of the honeycomb lattice H with a n x m fundamental domain. 

The fundamental domain of this periodic planar graph is represented on figure Q for n — m ^ 3. 
The vertices of the quotient Hi by the action of I? by translation are labeled by two integers, i and 
j ranging respectively from 1 to n, and from 1 to m. The weights of the edges around the white 
vertex labeled by (i, j) are denoted by aij, bij and c^ . By an appropriate gauge transformation 
one can spread the factors z and w in the magnetically altered Kasteleyn matrix K{z,w) so that 
the coefficients of this operator are Oy, hijW~^/^ and CijZ^^"^. The reference dimer configuration 
we will use is the configuration containing a-edges. 

One distinguishes three special classes of dual cycles in Gi , say A, B and C that cross only edges 
of a given type (respectively a, b and c). The lifts of these classes to H are represented on figure 
||. The B class is constituted by n vertical straight paths with homology class (0, 1), whereas the 
C class contains the m horizontal straight paths with homology class (1,0). The A class contains 
d = gcd(ri, m) paths with homology (^, ^). The three classes of cycles lift to H, forming three 
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(a) A-paths 



(b) B-paths 



(c) C-paths 



Figure 8: the three classes of possible frozen paths: A, B and C. 



classes of parallel families of straight lines. 

The Newton polygon of the weighted dimer graph G is a right triangle. Each side of the triangle 
corresponds to Gibbs measures for which all the paths of one of the three classes are frozen. The 
horizontal side contains n + 1 lattice points. The amoeba of the associated spectral curve exhibits 
n vertical tentacles separating n + 1 unbounded complementary components - solid phases in the 
phase diagram - each corresponding to a lattice point of the horizontal side. 

The frozen configurations are obtained as follows: when the Gibbs measure's slope lies, say, on 
the horizontal side of the Newton polygon, there is almost surely no edge crossing the B class of 
paths. These paths delimit thin strips where one sees either an infinite succession of a-edges or an 
infinite collection of c-edges. On G\ and in presence of a magnetic field (_B,i., By), the associated 
patterns have weights Hlii ^ij ^■^d e^" Hiii '^ij ■ The patterns with the highest weight correspond 
to the type of columns appearing in the configuration with probability 1. When the horizontal 
component of the magnetic field is very negative, the weight of the 'a' patterns are greater than 
the 'c' ones, but as increases, the weight of the second pattern becomes more important, and 
at some point, it becomes bigger that the first one. In the graph G the a-edges that were filling 
the space between the two B paths switch to c-edges. Generically the values of B^ corresponding 
to such a switch are all different. They correspond to the abscissae of the vertical tentacles of the 
amooeba. 

In a fixed window and when By is very largo, one sees columns of edges of the same type (a, or 
c) with a probability close to 1. When the magnetic field lies in a tentacle of the phase diagram, 
the system hesitates between two states for a given type of columns. With a probability p bounded 
away from and 1, the column is filled with o-edges, and with probability 1 — p, it is filled with 
c-edges. If one rcscalcs vertically the graph in the same time as By goes to +oC', then one will be 
able to see the transition between these two possibilities: between the two types of fillings, a &-edge 
is inserted. The edge creates a defect in the neighboring column that is supposed to be frozen. This 
discussion is quantified in the following proposition: 

Proposition 3. // 

m 

S.<^log^, (78) 
i=i ^^1° 
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then with probability going to 1 when By goes to infinity, the columns of type jo will be filled with 
a-edges. 

If the inequality is strict in the other direction, then they will be filled with c-edges with probability 
going to 1 when By goes to infinity. 

Proof. When the vertical component of the magnetic field is very negative, the main contribution 
to the characteristic polynomial is given by the configurations on Gi that contain no 6-cdgcs. One 
can choose to fill each strip of G\ between two consecutive frozen cycles either by a-edges or by 
c-edges. Choosing a filling with c edges induces a height change of m in the vertical direction. 
Therefore 

n / m m \ 

P^(z,w;)=P(e^-z,e^='u;) = n ma,, -(-l)"e^-znc,, +0(t) (79) 

j=i \i=i i=i / 

where t = e^" is small. 

Let b and w two vertices in the same strip jo in Gi . Denote by hy and wq lifts of b and w in 
the same column of G, separated by y fundamental domains. The entry of the inverse Kasteleyn 
operator between these two vertices is easily evaluated: recall that Qbw is the characteristic polyno- 
mial of the graph G where all the translated of b and w, as well as all the edges connected to these 
vertices, have been removed. Repeating the argument given above, we find the main contribution 
to it, 

n / m m \ 

Qi^{z,w) = Qy,^{e''-z,e''yw) = M^^^z^ H (H^ij ' (-^r^'"'' ^li^'ij] + 0{t) (80) 

j=i Vi=i i=i / 

where Afbw-z'^ is the weight of the dimer configuration of the strip jo of Gi deprived of b and w. The 
coefficient of the inverse Kasteleyn operator corresponding to these two vertices whose fundamental 
domains are separated by the lattice vector {x, y) is given by 

Kb (W, wo) - JJ^^ (81) 



Is 



z-y+^Mb^ dz dw 

't2 iUZi «yo - (-l)'"e^=^^ nlli Cijo) 2i7rz 2iTrw 
z-y+^Mb^ dz 



0{t) (82) 



+ 0(t) (83) 



/si (nlliOuo - (-l)'"e^"'^nlliCyo) 2i7r^; 
Suppose that B^ < Y]"^-, log the other case is similar. In that case, the pole located at 



nm 
i=l '^ijo 

nni 



is outside of the unit disk. If y — ^ < 0, then there is no pole at all in the unit disk, and therefore 
the integral over z is zero. However, i{y — d> 0, then the integral equals the opposite of the residue 
at zq, and 
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When bj^ and wq are the ends of an edge with weight aigj„ , Mbw = Ili^jo '^ijo ^^'^ U ^^nd 5 equal 
0. It follows that 

K-(b„wo) = =J^ + 0(t) = ^ + 0(t) (86) 

Thus, the probability of this a-edge, given by ai(jj(,K^^(bj,, Wq) goes to 1 when t goes to zero. 

On the other hand, in the case when b and w are the ends of a 'c'-type edge, then either y = —1 
or (5 — 1. In both cases, K^^(by, wq) is 0(t). Thus the probability of this edge goes to zero when 
t goes to zero. Such an edge is called non typical. □ 

A sequence of non-typical edges in a frozen column is initiated by the presence of a bead (a 
'6'-edge) crossing a neighboring wire. The analysis we made of the inverse Kasteleyn operator allows 
us to determine the distribution of the sequence of non typical edges in a frozen column. 

Proposition 4. The length of a succession of non typical edges in a frozen column has a geometric 
distribution. The parameter of the geometric distribution has an explicit expression in terms of 
ratios of lengths of dual edges. 

Proof. We suppose that in the frozen column jo, we only see a-edges with probability close to 1. 
The following inequality 

m 

B,<Y,\og^ (87) 

^=l ^^30 

is satisfied. Since we work only in one column, we will drop the index jo for the sake of simplicity. 
Given that the edge e = (w, b) with weight big is present in the dimer configuration, we compute the 
probability of seeing N successive c-edges after this bead. Denote by e^j; = (wi, bi), . . . , e^(wAr, bAr) 
the N c-edges. The weights of the edges around vertex are aj^] , , C[,j] , where [i] = (iq + i 
mod m) + 1. The conditional probability we want to compute is give by the following formula 



e G e and Vi = 1, . . . , TV e e e e 6 = 

P[e e e and Vi = 1,. ..,iV e e] 



' [e e e] 



^ N 



det Ajv+i 
K-i(b,w) 



where A^v+i is the following square matrix whose entries are inverse Kasteleyn operator coefficients: 





' K-i(b,w) 


... K-(b„w) ...1 










K-i(b,w,) 


K-(b„w,) 









(89) 



Since w and the white vertices stand on different sides of a frozen path, the associated coefficients 
K"^(b,Wi) are 0(t). More precisely, from (|65|), one has 



K-'(b,w,:) = tpbw. +0(t^ 



(90) 



Besides, if ? < j, the power of z in the numerator of ( pl| ) is positive, and it follows from computations 
made above that K"^(bj, Wi) is also 0(t). 
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For i g {1, . . . , N} and with the convention that wq = w, Wi_i and are the ends of an 
a-edge with weight a[i_i]. The same computations as above show that K^^(bi,Wi_i) = ^ ^ _^ ■ As 
a consequence, the asymptotic expansion of the determinant of Ajv+i is given by the product of 
these elements just above the diagonal times the last element of the first column. 

^ 1 

det An+1 ^ tpbw„ n + 0(t2) (91) 

As the probability of the bead is fcioPbwt + O(t^), the conditional probability we want is given by 

N 

P [iV successive c-edges |bead] = TT + 0(t) (92) 

Pbwjv 

Using proposition ^, one can rewrite p-h-w ^ I P\3wi as the following telescopic product 

N 

Plugging this into (p2|), one gets 

^ c U V ^ £(c ) 

P [N successive c-edges |bead] = TT "^'^ J' + 0(t) = TT , + 0(t) (94) 

T=~i "[^-ij'^bi '^-i T=i 

where £(a[i_i]) and i{c[i]) are respectively the length of the dual edges with weight a[i_i] and C[i] 
given by the mapping described in lemma |3| In particular, in the limit t ^ 0, the probability that 
the length L of this succession of non typical edges exceeds p fundamental domains equals 

Thus, in the limit, L has a geometric distribution. □ 




Figure 9: Illustration of proposition ^ Here is represented a frozen column of (horizontal) a-rhombi, perturbed by 
the presence of a defect (the 6-rhombus corresponding to the edge (w, b)), followed by a finite sequence of c-rhombi. 

Pushing further the above computations of the lengths of the non typical sequences of edges in 
frozen columns, one can derive the following: 
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Proposition 5. The limiting bead models on the different families of threads Bj are perfectly 
correlated: the distance between beads on each side of a frozen column converges in probability to 
zero. 

Proof. The details of the proof are omitted here, but by looking carefully at the determinants in 
the proof given above, one can in fact see that, for every t, the probability that a sequence of 
non typical edges in a frozen columns exceeds, say, is of order g^/^, and thus decays very fast 
when t goes to zero. Thus, in the vertically rescaled graph, the distance between two beads at the 
extremity of a sequence of non typical edges is close to 0. In the scaling limit, the distance between 
these beads converges in probability to 0. □ 



As a consequence, the picture of a typical dimer configuration for a Gibbs measure corresponding 
to a point in a tentacle of the phase diagram looks like the figure Eol 




Figure 10: A sketch of a typical tiling for a Gibbs measure corresponding to a point inside a tentacle. The black 
beads mark the transition between the two dominant type of edges. Next to them, beads also appear on the other 
threads to compensate the defect created in frozen columns. 



Acknoledgments 

This work has been initiated when the author was at the University Paris XI. The redaction 
was finished during a project at CWI, financially supported by the Netherlands Organisation for 
Scientific Research (NWO). 

We are grateful to thank Richard Kenyon for the many fruitful discussions and precious advices 
on this problem. We would like to thank also Daniel Slutsky and Jon Warren for pointing me to 
reference Q , as well as Scott Sheffield for showing me the random surface interpretation of the bead 
model. 



28 



References 

[1] J. Baik, T. Kriecherbauer, K. D. T.-R. McLaughlin, and P. D. Miller, Uniform 
asymptotics for polynomials orthogonal with respect to a general class of discrete weights and 
universality results for associated ensembles: announcement of results, Int. Math. Res. Not., 
(2003), pp. 821 858. 

[2] Y. Baryshnikov, GUEs and queues, Probab. Theory Related Fields, 119 (2001), pp. 256-274. 
[3] C. BOUTILLIER, Pattern densities in fluid dimer models, (2006). 

[4] D. J. Daley and D. Vere- Jones, An introduction to the theory of point processes. Springer 
Series in Statistics, Springer- Verlag, New York, 1988. 

[5] K. Johansson, Non-intersecting paths, random tilings and random matrices, Probab. Theory 
Related Fields, 123 (2002), pp. 225-280. 

[6] R. Kenyon, Local statistics of lattice dimers, Ann. Inst. H. Poincare Probab. Statist., 33 
(1997), pp. 591-618. 

[7] , The Laplacian and Dirac operators on critical planar graphs, Invent. Math., 150 (2002), 

pp. 409-439. 

[8] R. Kenyon and A. Okounkov, Planar dimers and Harnack curves, (2003). 

[9] R. Kenyon, A. Okounkov, and S. Sheffield, Dimers and amoebae, (2003). 

[10] M. L. Mehta, Random matrices, vol. 142 of Pure and Applied Mathematics (Amsterdam), 
Elsevier /Academic Press, Amsterdam, third ed., 2004. 

[11] G. MiKHALKiN, Amoebas of algebraic varieties and tropical geometry, in Different faces of 
geometry. Int. Math. Ser. (N. Y.), Kluwer/Plenum, New York, 2004, pp. 257-300. 

[12] G. MiKHALKiN AND H. RuLLGARD, Amocbas of maximal area, Internat. Math. Res. Notices, 
(2001), pp. 441-451. 

[13] M. Passare AND H. RuLLGARD, Amoebas, Monge-Ampere measures, and triangulations of 
the Newton polytope, Duke Math. J., 121 (2004), pp. 481-507. 

[14] M. Prahofer and H. Spohn, Current fluctuations for the totally asymmetric simple ex- 
clusion process, in In and out of equilibrium (Mambucaba, 2000), vol. 51 of Progr. Probab., 
Birkhauser Boston, Boston, MA, 2002, pp. 185-204. 

[15] S. Sheffield, Random Surfaces: Large Deviations Principles and Gradient Gibbs Measure 
Classifications , PhD thesis, Stanford University, 2004. 

[16] A. SOSHNIKOV, Determinantal random point fields, Uspekhi Mat. Nauk, 55 (2000), pp. 107- 
160. 



29 



